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Abstract 

We introduce a general expression which enables the parton distribution, unintegrated 
over the parton transverse momentum, to be obtained from the conventional parton den- 
sities. We use the formalism to study the effects of the transverse momentum qt of the in- 
coming partonic system on the calculation of the transverse momentum spectra of prompt 
photons produced in high energy pp and pp collisions. For the purposes of illustration, we 
use the double logarithm approximation. For large qt we calculate the effect directly from 
the perturbative formalism, whereas for small qt we bound the effect using two extreme 
hypotheses. In both qt domains we find that the shapes of the prompt photon spectra are 
not significantly modified, although the cross sections are enhanced. 



1 Introduction 

The cross sections for hard hadronic processes are conventionally described in terms of universal 
parton distributions a(x, /i^), with a = xq or xg, convoluted with the cross sections of the 
partonic subprocesses calculated in perturbative QCD at some large scale /i characteristic of 
the subprocess. For example, may be the transverse momentum of an outgoing parton or 
the mass of a heavy quark should it participate in the subprocess. In this paper we present 
the formalism which enables the parton distributions fa{x, k^, /i^), unintegrated over the parton 
transverse momentum kt, to be constructed from the conventional integrated distributions 



a{x, A^). We will see that it is important to distinguish the scale A from the hard scale fi. We 
first give the general prescription, and then we take the limits which correspond to the leading 
log(l/x) BFKL approach |I| and to the double log DDT prescription 0. 

These unintegrated parton densities fa{x, k'j^jj?) should be used to describe any hard had- 
ronic process, such as heavy quark production or the production of large Et jets 0. The only 
exception is deep inelastic scattering (DIS). If the kt integration is performed in the convolution 
of the unintegrated parton distributions with the deep inelastic partonic subprocess then the 
remaining x convolution is in terms of just the conventional parton distributions]^. 

To demonstrate the use of the formalism we study a particularly relevant and topical process, 
namely prompt photon hadroproduction. The production of prompt photons in high energy 
pp collisions has long played a key role in constraining the gluon distribution of the proton 
at large x through the dominance of the subprocess gq 7g. The reason is that the gluon 
enters at leading order, unlike its contribution to the description of deep inelastic scattering. 
However the description of the transverse momentum spectrum of the produced photons is more 
problematic than inclusive deep inelastic scattering. The observed pt-y spectrum in pp —>■ 7X 
(or pp — > ■yX) is found to be steeper than the prediction of perturbative QCD |Q, |^. The 
explanation of this discrepancy is usually attributed to the intrinsic transverse momenta kt 
of the incoming partons, which are usually assumed to have a Gaussian-like kt distribution 
IP, 1^. Thus part of the observed pt^y comes from the initial partonic kt such that the hard 
subprocess singularity da/dt ~ ^/Pt-y is approached more closely, and hence leads to a steeper 
Pty spectrum. However in order to describe the observed spectra it is necessary to introduce a 
kt spectrum with an average value which increases from (kt) ~ 0.5 GeV to more than 2 GeV 
as the collision energy ^/s increases from UA6, E706 p, ^ to Tevatron |]T0[ energies. Such large 
partonic (kt) cannot originate solely from the large distance confinement domain, but must also 
have a significant perturbative QCD component. Indeed it is easy to show that perturbative 
QCDQ gives a non-zero kf with a distribution in which (kt) increases as ^/s. Ref. |]T^ also 
contains a discussion of the partonic kt generated by perturbative QCD, but comparison with 
the pt^ data is based on a phenomenological Gaussian kt distribution. 

Here we study the role of both the non-perturbative and perturbative components of the 
incoming partonic kt in describing the observed pt.y spectrum in high energy prompt photon 
hadroproduction, that is pp (or pp) jX. Perturbative QCD is applied, at leading order, in 
the domain of /ar^fe partonic transverse momentum kt > qo, where go ~ 1 — 2 GeV is the starting 



point of DGLAP evolution in the global parton analyses [0, [13|. In the small kt < qo region 
it proves sufficient to consider two extreme hypotheses for the normalised non-perturbative or 
intrinsic kt spectrum: 



somewhat analogous example is the DIS renormalization scheme, where the higher order corrections are 
present in any process except DIS, which is used as the reference process. 

^Another correction to the perturbative form of the pt-f spectrum comes from the resummation of log(l — xt) 



terms, where xt = 2pty/y/s 11|. This effect only changes the spectrum at the larger values of xt, close to the 
kinematic boundary. 
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(a) a S{k^) distribution, i.e. no intrinsic kt from the confinement region, 

(b) a wide Gaussian distribution of width {k^) ~ (^q). 

From these extreme possibihties we conclude that the intrinsic partonic transverse momentum 
has a rather small effect on the pt^ spectrum. Moreover we find the effects of kf from the 
perturbative domain {kt > go) do not significantly change the pt-y spectrum. 

In Section 2 we introduce the unintegrated parton distributions and the parton-parton 
luminosity function Cab{xi,X2, qt), which is the probability of finding the incoming a, b partons 
(with momenta specified by Xi, kf^ and X2, k^t) with net transverse momentum = + kt2- 
In Section 3 we apply the formalism to the production of prompt photons in high energy pp 
or pp collisions. We convolute the luminosity functions with the hard gq jq and qq — > jg 
subprocess cross sections, taking into account the boost of the final 7+parton system coming 
from the incoming partonic qt- In Section 4, we show the effect of both the perturbative and 
non-perturbative components of qt on the predictions for the pt^ distributions at different ^/s, 
together with the measured pt-y spectra. Finally in Section 5 we present our conclusions. 



2 Unintegrated partons and the luminosity function 

To explore the effects of partonic transverse momenta we need to know the distribution of a 
parton a, say as a function of its transverse momentum kt, as well as of x. A straightforward 
way to obtain such a distribution is to consider the DGLAP evolution^ 



da as 



ainA2 27r 



(1) 



where the (integrated) parton density, a(a;, A^), denotes xg{x,\'^) or xq{x,\^). The first term 
on the right-hand-side describes the number of partons 5a emitted in the interval \^ < k'j < 
+ 5A^, while the second (virtual) term reflects the fact that the parton a disappears after 
the splitting. The second contribution may be resummed to give the survival probability Ta 
that the parton a with transverse momentum kt remains untouched in the evolution up to the 
factorization scale. The survival probability is given by the double logarithmic Sudakov factor 
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Uk,. ,) ^ e.p (- ^ f E . (2) 

Thus the probability to find a parton a with transverse momentum kt which initiates our hard 
process, with factorization scale fi, is 

fa{x,kl^^') = P„,,(^)a'g,fc2) dz) Ukt^fi). (3) 



■^For the g ^ gg splitting we need to insert a factor z' in the last integral of to account for the identity 
of the produced gluons. 



Now we have to specify the value of the infrared cut-off A, which is introduced to protect 
the 1/(1 — z') singularity in the splitting functions arising from soft gluon emission. In the 
original DGLAP equation (|I|) for integrated partons this singularity is cancelled between the 
real and virtual contributions. However after the resummation of the virtual terms the real 
soft gluon emission has to be accounted for explicitly as it changes the transverse momentum 
of the parton. That is we have to find the physically appropriate choice of the cut-off A, which 
comes from the coherence effect 
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The most convenient way is to go to the Breit frame 
where the angular ordering condition becomes the requirement that no soft gluons are emitted 
in the backward direction. In this frame the parton energy E = fi. Then the integral over 
the soft gluon momentum dz'/{l — z') = duj/uo covers the interval kt < uj < E, where u is the 
energy of the soft gluon. That is we take A = kt/ E = kt/ fi. Of course the same A must be 
used in both the real emission integral in and in the survival probability T in (^. Below 
we consider two limits of the general expression for the unintegrated parton density. 

In the leading log(l/x) or BFKL limit the virtual contribution in the DGLAP equation, 
(|l|), is neglected and the survival probability T = 1. Hence we come back to the familiar 
prescription 

2, d{a{x,\')) 



fa{,Xj kt 



(9 In A2 



(4) 



A2=fc? 



Note that in this limit the unintegrated parton density does not depend on /x, since there are 
no log(l/x) terms in the scale dependent part. 

In the double logarithm limit, @ may be written 



with survival probability 

Ta{X,fi) = exp 



d 



(9 In A2 



a{x, A^) Ta(A,/i) 



(5) 



\=kt 



A2 



c^siPt ) dp: 
27r 



Pt 



dz' 



-,2Ca 



(6) 



where Cq = Cp = {N^ — l)/2Nc and Cg = Ca = Nc- In this limit the cut-off A is numerically 
small so one can choose the value of A in (Q) equal to the value of A in (^). Therefore the 
second term of the right-hand-side of cancels the derivative dTa/ din X'^. To be precise, the 
only double logarithmic contribution comes from the singular 1/(1 — z') part of the diagonal 
splitting function Paa{z') in (|I]) which exactly equals the last integrand 2Ca/(l — z') in (^. 

To obtain the parton-parton luminosity function we have to perform a convolution over the 
distributions of the two incoming partons a, b 



Cabixi,X2,qt) = J faixi,kl,iJ^)fh{x2,kl,ii^)5^'^\ki + k2-qt 



d'^kid'^k2 



(7) 



where, for simplicity, we have omitted the t subscript on the ki. The luminosity element, Cabdqf, 
gives the probability that the incoming partons have a net square transverse momentum in the 
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interval {q^, + dq^). To leading order all the transverse momenta are strongly ordered, and 
so either ki <^ k2 — qt or ^ A;i — g*. In the first case we may integrate over ki giving 



Tiki 



a(xi, q^) Ta{qt, /i) fb{x2, /^^^ 



(8) 



and vice versa for ^2 -C ki. Thus the sum of both contributions may be expressed in the 
compact form 



Cab{xi,X2,qt) 



d 



(9) 



A similar expression for the Drell-Yan process was originally obtained in the classic work of 
ref. §. 



3 The prompt photon pt^ distribution 

To obtain the cross section for inclusive prompt photon hadroproduction we convolute the hard 
subprocess cross sections for gq •yq etc. with the parton luminosity of @. We obtain 

da _ da _ [ ^, dx,dx2 o 

-t^7 Is — ~ 1 TY' ~ / ^[Xi^X2, qt) aqt 

d'^p^ drj^irdpf^ J Xi X2 

X^di^ 5{V, - . . .) 5ipl -...)9{t,^- q^) e{\i\ - ql) (10) 

where is the azimuthal angle between and the photon transverse momentum p[ in the hard 
subprocess. That is p[ = |Msin^, where 9 is the polar scattering angle in the hard subprocess 
and = s is the invariant mass squared of the produced jq system. We thus have 

t = -— (1-cos^) (11) 

in the hard subprocess cross section da/di. The theta functions in ( [I0| ) impose the physical 
requirements that the partonic transverse momentum should be less than, first, the scale 
and, second, the momenta in the hard subprocess. The last condition comes from the fact 
that at leading order all the transverse momenta along the chain are strongly ordered. For 
the hard subprocess initiated by the gluon {ki) and the quark (^2) of Fig. |I] the virtualities 
k\^ k\ -C It means that q^ ~ max(fc^, A;|) < |£|. Fig. |l| also describes the qq 'jg subprocess 
but now /c| = i' plays the role of the hard momentum transfer and strong ordering implies 
<S |t'|. Thus the product Cda/di in ([T0|) should be understood as the sum over the 
different subprocesses, that is, as ^abdaab/diat where igg = i and i^q = i' in Fig. |I|. 

The introduction of the delta functions in (|TUp is simply a technical device to enable the 
differential cross section to be expressed in terms of r]^ and pf^. Let us study the relevant 
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Figure 1: A schematic diagram describing both the subprocesses gq 7g and qq '-/g, in 
which the hard momentum transfer squared is either t or t' respectively. 

kinematics. If the outgoing 7+ parton system, of mass M and rapidity rjM, is boosted with the 
incoming "partonic" transverse momentum qt, then the rapidity f]^ and transverse momentum 
Pt^ of the photon become 

1 f M^ + qs + Mcose \ 

Vy = W + o In T^— 77 a ' (12) 

2 \ M± + qs — M cos J 

. M!.„^,,M±«^^^iSi, (13) 

where 

qs = qt sine sin0, Ml = M^ + qt. (14) 

Equations ([T^ ) and (p!3|) specify the missing parts of the delta functions in (p^O]). As the variables 
rjM and M have been used to specify the Lorentz boost it is convenient to write the integration 
over the luminosity in ( |T0| ) in terms of dM^drjMdq^. This is most easily done by noting that 



'^(M^ - x^X2S + ql) 5 {vm - \ Hxjx^)) = (15) 



We may use the first delta function in ( |10|) to perform the drjM integration and the second delta 
function 5{pf^ ~ ■ ■ ■) to perform the integration. Then, using d\t\ = ^M'^ sin 9d9, we can 
rewrite ( p!0D in the form 

f r/n4'2 2\ ^It ^'P d^T 2M'^d9 n ( 2 2\ n ( \-\ 2^ 
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with M_L specified by 5{j)l —■■■), that is by 



sin (h 1 



= -qt—^ ± — JAp^ - ql cos^ 9 cos^ (17) 
Of course we may only use solutions of this latter equation which satisfy Mj_ > qt. 



4 Effects of partonic qt on the photon spectrum 

It is informative to look at the form of the luminosity function Cgq as a function of the transverse 
momentum qt of the incoming partonic {gq) system, before we consider the effects of qt on the pt-y 
spectrum of hadroproduced prompt photons. Fig. |^ shows Cgq obtained from (||), for values of qt 
in the perturbative domain {qt > qo), at the energies of the UA6 and CDF experiments. We use 



MRS(R2) partons |T5| and take q^ = 1.25 GeV^. At the higher Tevatron energy, ^/s = 1.8 TeV, 
we see that the qt distribution is extensive, whereas at the lower energy, ^/s = 24.3 GeV, the 
luminosity goes negative for sufficiently large qt- At this point we set £ = for larger values 
of qt- The negative values are an artefact of the double logarithmic approximation on which 
(|^) is based; the derivative dTa/dlnX^ does not exactly cancel the second (virtual) term in the 
DGLAP equation (|l]). The use of the full treatment based on (^ and (|^ would not suffer from 
this defect. 

We must, in addition, consider partonic transverse momentum in the domain qt < qo- In this 
region we cannot use perturbative QCD and, moreover, we have no knowledge of the parton 
distributions at scales A < go- On the other hand a significant part of the total luminosity 
comes from the region qt < qo for prompt photon hadropro duct ion at the lower energies, that is 
■\/s ~ 30 GeV. We may quantify this since the integrated luminosity Cq coming from the region 
qt < qo can be obtained directly from (|^) 

Co = a{xi,ql)Ta{qo,fi)b{x2,ql)Th{qo,iJ,). (18) 

Typical results are given in Table 1 in the relevant energy range. We show the luminosity 
integrated over the perturbative domain {qt > qo) together with Co for given values of xi = X2. 
We have fixed the value of xi = X2 = 0.3 in order to investigate the energy behaviour of the 
luminosity. For fixed Xi = X2 and increasing ^/s, the corresponding value of pt^ increases and 
we see that the fraction of the luminosity coming from the non-perturbative region decreases 
rapidly. At ^/s = 1800 GeV we also give the luminosity at a more experimentally representative 
value, Xi = X2 = 0.03, where it is much larger. 

In order to investigate the effect of the partonic qt on the prompt photon pt^ spectrum we 
need, not just the integrated luminosity Co, but also the qt distribution in the non-perturbative 
region. We take two extreme possibilities for the distribution in the region qt < qo- First, we 
assume that there is no intrinsic transverse momentum in the non-perturbative region, and set 

C{qt < qo) = Co 6{ql) (19) 
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Figure 2: The parton luminosity function Cgg{xi,X2,qt) of @ as a function of the incoming 
partonic qt for pt^ = 4 GeV at the UA6 energy of i/i = 24.3 GeV and for pt^ = 20 GeV at the 
Tevatron energy of ^/s = 1800 GeV. In each case we take Xi = X2 = xt = 2pt^/ ^/s. The hard 
scale /i is taken to be pt-y. 
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Table 1: The integrated luminosity from the perturbative {qt > qo) and non-perturbative 
(g* < go) domains for typical values of y/s, Xi, X2. We take q^ = 1.25 GeV^ and the hard scale 



(GeV) 


Xi = X2 


^gq{(lt > qo) 


^gqilt < qo) 


Fraction from 

qt < qo 


25 


0.3 


0.062 


0.160 


0.72 


50 


0.3 


0.095 


0.048 


0.34 


100 


0.3 


0.0851 


0.0103 


0.11 


500 


0.3 


0.0444 


0.0001 


2 X 10-^ 


1800 


0.3 


0.0291 


1.4 X 10^6 


5 X 10-^ 


1800 


0.03 


1.16 


0.0016 


1 X 10-3 



with Cq given by (|l^). Second, we assume a wide Gaussian distribution of partonic qt, 



C-iQt < qo) = 'Co 



(20) 



normalized to the region qt < go- From the results shown in Fig. ^ it turns out that the two 

for the non-perturbative qt distribution lead to a similar pt-y 



extreme hypotheses, (|19]) and 
spectrum. Insight into why this result may be anticipated can be obtained from the following 
simple estimate. 

On average the prompt photon gets one half of the qt of the incoming partonic system. Thus 
the effect of the partonic qt on the hard subprocess cross section is 



da 
di 



1 



Ptj - 2lt 



4 

Pt7 



1 + 



3(Pt7 • Qt) 
4 

Pt7 



'2p% 



4 

P*7 



1 + 4 

Ptl 



(21) 



after the angular integration. So the relative difference between the qt < qo effects calculated 
using (|1^) and (pO|) is expected to be only about gQ/2pj^, since in the non-perturbative domain 

{ql) < ql 

At this point one may ask why the traditional phenomenological treatments |^, |^ of the 
intrinsic partonic qt have found such a large effect. The reason is that the full kinematics were 
not taken into account. Instead of giving the Lorentz boost to the produced 7+parton system, 
the photon on its own received the total partonic qt- To reproduce the same effect with the 
correct kinematics would require about twice as much partonic qt- Moreover strong-ordering 
of transverse momenta should be imposed which prevents the hard subprocess approaching the 
t ~ region. 
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Figure 3: The effect of partonic qt arising from tlie non-perturbative {qt < qo), perturbative 
{Qt > Qo), and total contributions to tlie pt-./ spectrum of prompt plioton production in pp 
collisions at ^/s = 24.3 GeV and in pp collisions at 1800 GeV. Two extreme hypotheses are 
used for the non-perturbative contribution, but the predictions are very similar (and in fact are 
indistinguishable on the ^/s — 1800 GeV plot). We see that the non-perturbative contribution 
dominates at ^/s — 24.3 GeV, whereas the perturbative dominates at ^/s — 1800 GeV. The 
hard scale /i is taken to be pt-y and Qq — 1.25 GeV^. 
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It is informative to discuss the results at the lowest (UA6) and the highest (Tevatron) 
energies separately. At the lowest energy we see from Fig. |^ that the contribution from the 
perturbative {qt > go) region is relatively small. Thus, even though it depends on the choice of 
the hard scale /i, the ambiguity is not so important. Moreover we note that the two extreme 



choices of the non-perturbative qt distribution, (p!9D and (|20|), give only a 20% uncertainty in 
the pt-y spectrum. That is partonic qt smearing has little effect on the pt^ spectrum, as may 
be anticipated from (|2T|) . However there are other related points to consider. We see that 
the formalism requires that the parton densities are sampled at scale go, which enhances the 
cross section in comparison with the conventional leading order treatment in which the partons 
are sampled at the hard scale /i. The qt < qo contribution depends on the scale fi through the 



survival probabilities Tq and Tb in Cq of (|T^) . If, for illustration, we take the hard scale = cpt^ 



then, to double log accuracy with fixed as, we have 

Ta ^ exp ^ asicpt,) (^^^^ . (22) 

Thus an increase in scale, that is a larger c, makes da/dpt^ steeper and the cross section smaller. 
The prediction for the prompt photon spectrum, da/dpt.y, at the Tevatron energy is much more 
certain. Here the contribution for the non-perturbative qt < qo domain is negligible (see Fig. 
and so the calculation is more under control and is less scale dependent. The difference between 
using (H) and using an "unsmeared" luminosity 

Cfe = S{qt ) a{xi, fi"^) b{x2, fi"^) (23) 

is not great. The difference is due to the qt smearing itself and also due to the fact that in (|^) 
the partons are sampled at qt. This latter effect was not incorporated in the existing treatments 
of partonic smearing of the pt^y distribution. 

To demonstrate the sensitivity to the choice of the hard scale, in Fig. ^ we compare the 
results for fi = Pf^/2 with those for fi = pt^ for three typical energies at which experimental 
measurements exist. A third choice of the hard scale, n = M/2, gives similar results to those 
obtained with /x = pt^. For comparison the unsmeared cross sections are shown by dashed 
curves. Very little of the difference between the continuous {qt ^ 0) and the dashed {qt = 0) 
curves is directly attributable to the introduction of the partonic transverse momentum qt. As 
mentioned above, the main difference is due to the proper inclusion of the parton densities, 
which are sampled at scale qt (with qt > qo) in contrast to the conventional unsmeared result 
in which they are sampled at the hard scale /i. Of course in a realistic fit to the data we 
would need to implement the more precise form (|^) of the unintegrated parton densities, as 
well as computing the next-to-leading order expression for the hard subprocess cross sections. 
However for the largest energy, =1.8 TeV, where the pt-y is large and leading order is a 
good approximation, the theoretical uncertainties are seen to be small0. In fact the /i = pt.y 

*The difference between the predictions for /i = pty and /i = pt-y/i may be considered as an estimate of tlie 
size of the next-to- leading order contribution. 
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prediction gives a reasonable description of the CDF data, particularly at the larger values of 
Pfy, see Fig. |^. Figs. ^ and ^ also show that, in common with other analyses^ there is a large 
disagreement between the predictions and the pt^ shape of the E706 data. 

5 Conclusions 

This paper should not be regarded as an attempt to better fit data, such as the pt-y spectrum 
of prompt photon hadroproduction, but rather to illustrate the appropriate formalism that 
should be used to incorporate the transverse momenta of the incoming partons. The crucial 
point is that, in general, according to the /cj-factorization prescription 0, cross sections should 
be calculated in terms of parton distributions, fa{x,kf,fi'^), unintegrated over the transverse 
momentum. Such distributions can be constructed from the conventional parton densities, 
a{x, A^), as described in Section 2. 

Prompt photon hadroproduction is an ideal example to illustrate the necessity of using un- 
integrated parton distributions. For simplicity we chose to work to double logarithm accuracy, 
rather than the full treatment presented in Section 2. Contrary to the existing phenomenolog- 
ical treatments of the intrinsic transverse momenta of the partons, we find that the qt effects 
give themselves only a small modification to the da/dpt^ spectrum. Of course there could have 
been a large effect if qt ~ pt^^, which would have allowed the singularity of the hard subprocess 
amplitude to be approached. However the strong ordering of transverse momenta must be 
incorporated in the formalism, via 9{\i\ — q'f), which prevents this happening. In addition, our 
study reveals that the qt of the incoming partonic system should be correlated with the scale at 
which the parton densities are sampled. This fact affects both the shape and the normalization 
of the predictions, especially at the lower energies. 

An interesting result of our study of the parton luminosity is that, at low energy where 
the contribution from the qt > qo region is small, the factorization procedure is more direct 
and physically transparent. Recall that we have to introduce a factorization scale to separate 
the hard subprocess interaction (which is calculated perturbatively) from the universal parton 
distributions which originate from large distances (that is, from the confinement region). These 
incoming partons starting at the input scale go participate in the perturbative evolution which 
changes the original distributions. However for larger x (that is for the lower energies) we sample 
the parton densities directly at the low input scale qo, see Fig. |^ and (p!8|). From the form of the 
luminosity we see that essentially no new partons appear in the evolution to larger values of 
qt- The only perturbative effect is the introduction of the survival probability, T^, which has a 
simple physical interpretation and does not alter the kinematics of the partons sampled. This 
means that the low energy data|^ measure the input distributions directly. On the other hand 
at the smaller valuesQ of x the perturbative contribution qt > go dominates. In this case the 

^For example, see the (unsmeared) results in Refs. ||, ||, 0]- 
^For example, UA6 data at pt^ — 5 GeV. 
^For example, Tevatron data at pty = 40 GeV. 
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Figure 4: The scale dependence of the predictions for production of prompt photons in pp 
colhsions at a/s = 24.3 GeV, pBe colhsions at y/s = 38.8 GeV and pp colhsions at y/s = 1.8 TeV 
shown together with UA6 0], E706 [Q and CDF data respectively. The continuous curves 
are the predictions with the incoming partonic transverse momentum qt included, whereas the 
dashed curves correspond to the unsmeared {qt = 0) results in which the integrated partons are 
sampled at the hard scale /i. In each case the upper curve corresponds to the scale = Pt-yf^, 
while the lower corresponds to yU = pfy. 
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Figure 5: A clearer comparison of the data also plotted in Fig. ^ with our theoretical calculations 
at the scale /i = pt-y] we use xt = Ipt^j^ as the transverse momentum variable. 
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partons evolve from the input scale qo and are finally sampled at the hard factorization scale 
H ~ pt^. The cvohition of the integrated parton densities in the conventional approach modifies 
the X behaviour of the distributions assuming that the transverse momenta kf of the partons is 
small. The formalism that is presented here embodies the modification of the kt distributions 
of the partons during the evolution, as well as the modification of their x dependence. That 
is we are able to determine, perturbatively, the qt = + ^2*1 distribution of the incoming 
partonic system, which leaves no room for additional phenomenological qt smearing. 

In summary, we have presented the formalism to construct the parton distributions un- 
integrated over the parton transverse momenta, fa{x,k^, /j,"^), from the known conventional 
parton distributions. These distributions should be used to calculate cross sections initiated by 
hadrons. We applied the formalism to prompt photon hadroproduction in the double logarithm 
approximation. In this way we are able to quantify the effects of the partonic transverse mo- 
mentum, qt, on the transverse momentum, pt^, spectrum of the produced photons. We find that 
the correct treatment requires that the integrated partons are sampled at the scale qt. Clearly 
this is only meaningful in the perturbative domain qt > qo, but we are able to show that the 
uncertainty due to the contribution from the non-perturbative domain is only of relative size 
of the order of qQ/2p^^. 
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